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Abstract SuperCDMS is the next phase of the Cryogenic Dark Matter Search 
experiment, which measures both phonon and charge signals generated by par- 
ticle recoils within a germanium target mass. Charge signals are employed both 
in the definition of a fiducial volume and in the rejection of electron recoil back- 
ground events. Alternatively, phonons generated by the charge carriers can also 
be used for the same two goals. This paper describes preliminary efforts to ob- 
serve and quantify these contributions to the phonon signal and then use them to 
reject background events. A simple analysis using only one pulse shape parameter 
shows bulk electron recoil vs. bulk nuclear recoil discrimination to the level of 
1:10-^ (hmited by the statistics of the data), with httle degradation in discrimina- 
tion abihty down to at least 7 keV recoil energy. Such phonon-only discrimination 
can provide a useful cross-check to the standard discrimination methods, and it 
also points towards the potential of a device optimized for a phonon-only mea- 
surement. 

PACS numbers: 29.40.Wk 63.20.kd 



1 Phonons in a SuperCDMS Detector 

SuperCDMS (Cryogenic Dark Matter Search) Soudan is in preparation for instal- 
lation at the Soudan Underground Laboratory in northern Minnesota. The target 
mass consists of ^15 large roughly cylindrical (3.81 cm radius, 2.54 cm thickness, 
^600 g) ultra-pure germanium crystals. These target masses are instrumented 
on the flat surfaces with interleaved arrays of phonon- sensing Quasiparticle-trap- 
assisted Electrothermal- feedback Transition-edge- sensors (QETs) and charge- sensing 
electrodes. ^ ^ The measurement of both the phonon energy and charge energy of 
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each event enables the ratio of these two energies to be used as a powerful event- 
by-event discriminator between low-charge-production nuclear recoils (a possible 
dark matter signature) and high-charge-production electron recoils (the vast ma- 
jority of background events). 

The interleaving of biased electrodes (typically +2 V on one side, -2 V on the 
opposite side) with QET arrays (essentially at V) defines a ^0.5 V/cm drift field 
in the bulk of the target mass. Within ^1 mm of the top and bottom surfaces, 
however, the field created by the interleaved structures is both much stronger than 
the bulk (^20 V/cm), and also largely parallel to the surface. This field shape 
prevents carriers created in this region from propagating to the opposite side elec- 
trodes, and the resulting side-asymmetric charge signal tags such events as near- 
surface events'^ (near-surface electron recoils would otherwise form a dangerous 
background's). 

Much of the essential information content of the charge measurement is du- 
plicated in the phonon measurement, through the production of secondary phonon 
populations as charges propagate through the crystal and then enter aluminum 
structures at the surface. These different phonon populations are produced with 
different initial energy distributions. In Germanium, high energy phonons exhibit 
high isotopic scattering rates (t^^ = [36.7 x 10""^^] for Ge)^meaning that they 
have short mean free paths and propagate diffusively (i = [5.4 km /s]Ti = [1 .5 x 
10"^"^ m//] at long wavelengths). Only 6.1% of the top and bottom crystal sur- 
faces is covered with phonon-absorbing aluminum (largely in the form of QETs), 
meaning that, per surface interaction, a phonon is far more likely to reflect off a 
polished Ge surface than be absorbed in a QET. A phonon's timescale of surface 
absorption then is roughly proportional to its rate of surface interactions, and is 
therefor highly dependent on mean free path. Phonons with a mean free path on 
order of the detector size (ie, which propagate ballistically) are absorbed with the 
maximally slow decay constant of 755 fis, and higher energy phonons are ab- 
sorbed with significantly faster decay constants. Phonons anharmonically decay 
in energy (t~^ = [1.61 x 10~^^]v^)^, meaning that mean free path increases with 
time. This energy-dependent diffusion and decay process is typically termed "qua- 
sidiffusion". Here we summarize the relevant characteristics of the major phonon 
populations in a SuperCDMS detector, emphasizing where they begin in their qua- 
sidiffusive evolution: 

- Primary phonons produced by the recoil event itself are initially highly en- 
ergetic (v >1 THz, ie ^ < ^ lOO^im). If the event occurs far from a sensor 
surface, the diffusive behavior slows the arrival at the surface and lengthens 
the eventual mean free path, slowing the absorption rate at that surface. If 
the event occurs near a sensor surface, there is httle delay in arrival, and the 
comparatively short mean free path increases a the rate of absorption at that 
surface. 

- Neganov-Luke phononsP^ are created as charge carriers are drifted by the 
electric field. Such phonons are created very soon after the event time (the 
charge drift time is ^ 1 ;Us). In the low field of the detector bulk, Neganov- 
Luke phonons are emitted with low energies and correspondingly ballistic free 
paths, whereas the majority of Neganov-Luke phonons are produced in the 
strong near-surface fields at higher frequencies (^300 GHz and ^700 GHz 
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for electrons and holes) with correspondingly short mean free paths (^1.8 cm 
and ^600 iimj^^md fast absorption rates. 
- Relaxation phonons are emitted when charge carriers cross from the Ge into 
the superconducting Al surface structures, converting their Ge gap energy 
(0.75 eV) into Al quasiparticles (of energy 2Aai = 200 ^ueV) and a popu- 
lation of initially high energy phonons. In the presence of superconducti ng Al, 
such phonons experience rapid downconversion (at time scales of ^ 1 jisj^^^, 
meaning that relaxation phonons are ballistic at the time scales considered 
here. 

These three phonon populations combine together to produce the observed 
phonon pulses and their associated rising and falhng timing characteristics. We 
will from now on ignore relaxation phonons, since they are initially ballistic (their 
absorption time constant is a slow ^755 fis, and a minimal contribution to pulse 
shape). A recoil occurring in the crystal bulk will produce a roughly side- symmetric 
signal of near-surface Neganov-Luke phonons, which will arrive at both the top 
and bottom surfaces at fast times (^l < t < ^5 fis) and are absorbed rapidly by 
these surfaces. Primary phonons will reach the top and bottom surfaces more grad- 
ually, as quasi-diffusive propagation allows ( ^1 < t <^15 fis). These primary 
phonons will have shorter mean free paths (and will be more quickly absorbed) at 
the closer surface, and will have longer mean free paths (and will be more slowly 
absorbed) at the further surface. When an event occurs within the near- surface 
electric field region, it will be distinguished by a boosted population of Neganov- 
Luke phonons at one surface (with fast absorption times) and a lack of any distinct 
Neganov-Luke sharp rising edge in the opposite side pulse. 



2 Quantifying Phonon Pulse Characteristics in SuperCDMS Detectors 

The understanding of phonon populations given in the previous section predicts 
certain shape characteristics as a function of both proximity to the absorbing sur- 
face and quantity of charge carriers (and resulting Neganov-Luke production). 
Looking at recent calibration data (Figure 1), we do see such behaviors, confirm- 
ing our understanding. 

To make use of these phonon pulse shape variations as event type discrimina- 
tors(distinguishing electron recoil and near-surface recoils from nuclear recoils), 
we first reduce the raw traces to a small number of simple quantities, and then 
combine these quantities in a way that accommodates the large position-dependent 
variation. First, each of the raw traces from the detector's eight QET channels were 
fit using the following functional form: 

where the total pulse is treated as the sum of a "fast" pulse and a "slow" pulse, 
with rising and falhng time constants Rfast'> ^siow^ ^fast^ and Fsiow^ and start time 
offsets Tfast and Tsiow The fast pulse can be thought of as the contribution of 
highly-diffusive quickly-absorbing phonons, while the slow signal can be thought 
of as the contribution of more ballistic slowly-absorbing phonons, but the quasi- 
diffusive nature means that in reahty there is not a simple two-category distinction. 
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Fig. 1 Example phonon pulses from calibration runs with a SuperCDMS Soudan detector. The 
four side 1 QET channels are summed and shown in blue (dark grey); the side 2 channels are 
summed and shown in red (light grey). All eight events shown are of similar total phonon energy 
(~50 keV), and are located at different positions along the central axis of the cylindrical crystal 
(as determined by relative weighting among the eight phonon channels). Surface 1 and surface 2 
events were selected using an additional a requirement that the charge signals from the two sides 
be strongly asymmetric. 



The falling time constant Fsiow of the slow pulse was observed to be nearly iden- 
tical for every event, and was set to 755 ^us, the observed rate of absorption of the 
late-time uniform bath of low energy phonons. Each event's resulting 7 fit pa- 
rameters for each of 8 phonon channels were further reduced by summing the fits 
for each side and then finding key points (10%, 20%, 30%, etc.) along the rising 
and falling edges of these side- summed fits. 

The partition of energy between the eight channels also contains significant 
discrimination information (in addition to position information). The amplitude 
of each channel's pulse was obtained using an optimal filter, in which a template 
pulse of fixed shape was scaled to best match the amplitude of the (variable shape) 
phonon pulse. These measured amplitudes, then, are shape-dependent. Distribu- 
tions of some example partition and pulse shape quantities are shown in Figure 2. 



3 Discrimination Based on Phonon Pulse Quantities 

As a very first look at discrimination based on phonon pulse shape characteristics, 
we here show a simple example of bulk electron recoil (ER) vs bulk nuclear recoil 
(NR) discrimination using only one pulse shape characteristic: the 40%-to-70% 
risetimes. This simple discrimination example serves as a lower bound on the 
abilities one might expect when using a combined analysis of all phonon pulse 
quantities. 

We first construct ER and NR populations. Identical charge-based fiducial vol- 
ume cuts are enforced on both populations, and measured [charge: phonon] energy 
ratios are used to categorize events as either ER or NR. After these two popula- 
tions have been defined, a new pulse shape discrimination quantity is constructed 
using the 40%-to-70% risetimes (for side 1 and side 2) for each event: "radius" = 
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Fig. 2 Example phonon discrimination quantities. On the left, an example of partitioning be- 
tween channels is shown. The vertical axis is an axial partition quantity (sidel - side2)/(sidel 
+ side2), and the x axis is a similarly defined radial partition quantity. On the right, the time 
difference (in ^s) between the 40% and 70% points on the rising edge of the side 1 summed 
pulse (x axis) and side 2 summed pulse (y axis) are shown. In both plots, calibration events 
are colored using a charge-signal-based categorization as either bulk electron recoils (blue or 
dark grey), bulk nuclear recoils (green or light grey), or side 1 electron recoils (black). The total 
phonon energy for these events are between 7 and 20 keV recoil energy, where recoil energy has 
been scaled from total phonon energy using a nuclear recoil assumption for all events. 



([40%-to-70% side 1]^ + [40%-to-70% side 2]2)V2, One can see in Figure 2 that 
this combined quantity is largely position-independent. 

Figure 3 shows a histogram of this timing quantity, plotting only the lowest 
energies inspected in this analysis (7 < E^-ecoU < 20 keV, defined by assuming that 
all events are nuclear recoils and scahng the total phonon energy accordingly). Al- 
though the ER and NR distributions overlap somewhat, there are no slow ER out- 
liers to the statistics available in the calibration dataset. The right panel of Figure 

3 shows the NR acceptance fraction vs ER leakage fraction as one varies a ID cut 
threshold. Discrimination better than 1:10-^ is seen, and it is further seen that this 
discrimination shows no degradation with energy down to at least 7 keV Ey^QQU. 

4 Conclusions 

Although significant pulse shape discrimination ability has been demonstrated, 
phonon-only discrimination is only beginning to show its promise. The simple 
analysis performed here can be improved in many ways, and looking further on the 
horizon, this rough analysis using a detector optimized for charge-electrode-based 
discrimination suggests that a detector optimized specifically for phonon-only dis- 
crimination could take advantage of the extreme sensitivity of transition edge sen- 
sors to extend event-by-event discrimination well below the typical CDMS low 
energy threshold of ^10 keV (set by charge readout noise). Such phonon-only dis- 
crimination capabilities should be possible if detectors are specifically optimized 
for this goal, through the reduction of TES internal thermal fluctuation noise 
and also through an increase in pulse shape differences themselves by increasing 
the total phonon-absorbing Al area of the QETs. 
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Fig. 3 A simple look at pulse shape discrimination between bulk nuclear recoils and bulk elec- 
tron recoils, using 40%-to-70% risetimes. Nuclear recoils are green (light grey); electron recoils 
are blue (dark grey). For discussion, see text. 
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